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Abstract: Spontaneous loss of charge by charged black holes by means of pair-creation of charged 
Dirac particles is considered. We provide three examples of exact calculations for the spontaneous 
discharge process for 4D charged black holes by considering the process on three special non- 
rotating de Sitter black hole backgrounds, which allow to bring back the problem to a Kaluza-Klein 
reduction. Both the zeta-function approach and the transmission coefficient approach are taken into 
account. A comparison between the two methods is also provided, as well as a comparison with 
WKB results. In the case of non-zero temperature of the geometric background, we also discuss 
thermal effects on the discharge process. 
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1^. Conclusions 



1. Introduction 

Spontaneous loss of charge by a charged black hole is a relevant topic in the framework of quantum 
effects in the field of a black hole || . It belongs to the framework of phenomena which are due 
to vacuum instability in presence of an external field, with consequent pair creation. Quantum- 
electrodynamics effects in presence of an external electric field have been in particular a key-topic 
which has been extensively discussed. Being our interest oriented toward an application to black 
hole physics, we limit ourselves to quote two seminal papers |3[ ^ and, in the recent literature, 
Refs. |6[ 0|. An effective description of pair creation phenomenon for static charged black holes 
was provided by Damour, Deruelle and Ruffini in a series of papers [|[ ||, |l^^. To sum up, on these 
backgrounds the Hamilton- Jacobi equations (H-J) for a classical charged particle can be easily re- 
duced to quadrature by means of variables separation. In particular, the radial equation describes 
a one dimensional motion of a particle in a given potential. The H-J equation, beyond a positive 
energy potential, determines a negative energy potential which at classical level must be discarded. 
However, at quantum level, negative energy states must be included, and a quantum interpretation 
to this couple of potentials can be given. The positive energy potential determines the allowed pos- 
itive energy states, whereas the negative energy potential determines the allowed negative energy 
states. The usual separation of these states occurring in absence of external fields is not ensured a 
priori, and there can be regions where an overlap of positive and negative states for the particle is 

^Also rotating black holes were treated therein. 
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allowed, i.e. the Klein paradox takes place. In these level-crossing regions, by means of tunneling 
between negative and positive states, pair production of charged particles can take place with a rate 
determined by the transmission probability for the particle to cross the forbidden region between 
the two potentials, and can be computed e.g. in the WKB approximation. 

We improved this semiclassical picture in the case of anti de Sitter Reissner-Nordstrom black holes 
showing that the potentials have a direct interpretation at the quantum level without referring to 
the classical H-J equation Then, for the class of de Sitter Reissner-Nordstrom black holes we 
found that level-crossing is always present, due to the peculiar occurrence of both a black hole event 
horizon and a cosmological event horizon and we also considered a particular limit case, when 
the black hole horizon radius r+ equates the cosmological horizon radius rc- the Nariai black hole 
p3| , p^ , p^ . The aforementioned class of solutions contains further limit cases, corresponding to 
the extremal cases r_ — r+ = which are called ultracold solutions of type I and II ||l^, A 
careful WKB analysis was also performed for the Nariai case and the ultracold ones. 
Herein, we develop our analysis of the pair-creation process associated with the black hole electro- 
static field, and fully exploit the fact that the aforementioned special backgrounds allow an exact 
calculation of the vacuum instability. As a consequence, we can provide for the first time, to our 
knowledge, exact results for the instability of 4D charged black holes. We point out that our back- 
grounds are of a special character: in all the cases the geometry involved is the one of a Cartesian 
product Af^'^ X S'^ where M^'^ is a two dimensional spacetime and is a sphere with constant 
radius (there is not any non-constant warping factor). Moreover, the fluxes do not involve the 
sphere directions so that the sphere could be effectively considered as an internal space and the 
problem can be reduced to a two dimensional effective theory by means of a Kaluza-Klein reduction 
p6| , p^ , p9[ . Here we will recall manifestly the strategy of a K-K reduction in part of our analysis. 
Indeed, for all cases, the Dirac equation will be reduced to a two dimensional Dirac equation on 
the M^'^ background, with the two dimensional spinors obtained by the two dimensional reduction 
of the four components Dirac spinors and the mass spectra corrected by the K-K modes. We also 
recall that in a K-K reduction the latter corrections are provided by the harmonic analysis of the 
internal space S^: 

2 2 2 , \ 

where is the particle mass and A; are the eigenvalues of the Laplacian operator — A52 on the 
internal space. In our case we are involved with harmonic analysis for spinors, and then the Dirac 
"angular momentum" operator eigenvalues k = ±(j + ^) G Z — {0} appear in A;. 

We also take into account the fact that, both in the Nariai case and in the ultracold I one, 
the real quantum state to be considered is not the Boulware-like state corresponding to standard 
quantum vacuum, but the Hartle-Hawking state associated with the black hole temperature. Then, 
we discuss also how pair-creation due to the electrostatic field of the black hole superimposes to the 
thermal radiation effect which is present in the given backgrounds. We show that, in the thermal 
means of the Dirac field number operators, a standard thermal contribution appears together with a 
term which is still related to the aforementioned vacuum instability, except for a further dependence 
on the background temperature which is induced by the Hawking effect. 

The plan of the paper is the following. In section ^ we recall the transmission coefficient 
approach to vacuum instability, and consider also how the instability affects thermal states in the 
framework of Thermofield Dynamics. In section |^ we discuss the ultracold II case, both in the 
transmission coefficient approach and in the zeta function approach. In section |^, an analogous 
analysis is carried out for the ultracold I case. In section ^ the Nariai case is considered. In this 
case, the zeta function approach requires a recently developed calculation strategy |2^, which is 
sketched in the present paper. For all the cases a comparison with WKB results is done. In section 
1^ conclusions appear. 



- 2 - 



2. Vacuum instability and thermal states 



In this section, we first recall some aspects of the pair creation due to vacuum instability, in 
particular we focus on the so-called transmission coefRcient approach to the evaluation of the 
instability, which is associated with the presence of an imaginary part of the effective action 
which in the following is approached also by means of ^-function techniques. Then wc consider how 
instability in external fields affects thermal states, referring to Thermofield Dynamics approach for 
a general fermionic case, although our interest is in the black hole case. 

2.1 Vacuum instability in the transmission coefflcient approach 

There are several ways one can deal with vacuum instability. One consists in Schwinger's approach 
iQ, |2l|, with calculations carried out in the correct space-time signature, or in its ^-function 
variant where calculations are developed in Euclidean signature and then a rotation to real time is 
performed. As to vacuum instability, we adopt the ^-function variant, and a double check of our 
results is also provided by the so-called transmission coefficient approach. We recall shortly some 
aspects of the latter, in which, following H, it is also possible to reconstruct the probability 

of persistence of the vacuum. Let us introduce, for a diagonal scattering process, 

r,IN _ „OUT j,rp OUT 

where rii stays for a negative energy mode and pi for a positive energy one. Ti is the transmission 
coefRcient and Ri is the reflection one. Moreover, in one defines 

V^ ■■= (2.2) 

which coincides with the mean number per unit time and unit volume of created particles [|| . Cf. 
also 11^. In the case of fermions, the result is 

\R,\^^l-r,,, (2.3) 

which excludes any superradiant phenomenon for Dirac particles. As thoroughfuUy discussed in 
p^ , it is crucial to consider group velocity for the asymptotic behavior of the solutions, in order to 
get a correct physical result. In fact, an erroneous consideration of phase velocity in place of group 
velocity could easily lead to conclude that superradiance exists also for the Dirac case. 
By interpreting a la Stueckelberg the scattering process, one can also obtain 

n?^^ = R^'nl^ - Ri'T^p?^^, (2.4) 

which is interpreted as the scattering of a negative mode incident from the future and which is in 
part refracted in the past and in part reflected in the future. The new reflection amplitude —R~^Ti 
is such that the reflection coefficient 

\R7'n' = (2.5) 

can be interpreted as the relative probability for the creation of the pair nf^^.pf^'^. The absolute 
probability is obtained by multiplying the relative one times the probability pi,^ to get zero pairs 
in the channel i, and then the probability „ of n pair for fermions is 

P.,„=P.0(^3^. (2.6) 

The normalization condition 

1 

^Pz.n = l, (2.7) 
n=0 
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for fermions leads to 

P»,o = 1 - (2.8) 

The persistence of the vacuum is given by 

Po = Y[p^,o^exp{~2lmW), (2.9) 

i 

and then 

oo ^ 

21mW = ~Y.^og{l~r^.)^J2J2p^ (2.10) 

i i k=l 

for fermions. For bosons see e.g. 
2.2 Finite temperature effects 

In the case of ultracold I and Nariai manifolds, one has to take into account that we deal with a 
black hole manifold endowed with a non-zero temperature. As a consequence, we have to consider 
quantum instability not simply for a vacuum state which corresponds to the Boulware-like state 
of standard Schwarzschild solution, but for the state which plays the role of Hartle-Hawking state 
for the given solution. Since the discover of the Hawking effect, a very fine construction of the 
thermal state living on a finite temperature black hole manifold characterized by a bifurcate Killing 
horizon was introduced by Israel on the grounds of Unruh analysis pTf and of the thermofield 
approach to thermal physics introduced by Takahashi and Umezawa |28| , p9| , |30| . Israel discov- 
ered that the HH state (hence called also Hartle-Hawking-Israel state) corresponds to the thermal 
vacuum of thermofield approach with the temperature equal to the black-hole temperature, and 
that the would-be fictitious states of the thermofield approach correspond to the states in the left 
wedge of the extended solution (if one is living in the right wedge). In our case, one possibility 
is to consider an analogous construction; alternatively, we adopt a more "liberal" attitude, in the 
sense that we appeal to Thermofield Dynamics formalism and describe the equilibrium state (KMS 
state) without caring about the reality of the "would-be fictitious" states in a "specular wedge" of 
the extended manifold. One could also appeal to the approach developed in (where a KMS 
quantization, leading to a KMS state, is introduced, without any doubling of the Hilbert space). 
We show that the transmission coefficient approach at finite temperature is still a valid approach 
to analyze the quantum instability problem at hand. Thermofield dynamics also helps a straight- 
forward generalization of quantum instability to the case where the initial (in) and the final (out) 
states are thermal states (at the same temperature) instead than vacuum ones. We could as well 
start from results given in [^2| |3^ and also in ||2^, which analyze stability topics in quantum 
electrodynamics. In order to check if there is instability in the thermal state at the Hawking tem- 
perature, we adopt the following "standard" strategy: we calculate the mean number in the "in" 
thermal state of "out" particles in the ^-mode, and see if there exists any deviation from a purely 
thermal distribution. Equivalently, we could evaluate the thermal mean of the number of "out" 
particles in the Z-mode minus the number of "in" particles in the Z-mode (cf. |^) and then see 
the net efii'ect of a possible quantum instability. In what follows, our focus is to the case of black 
hole backgrounds with a single temperature, and then (3 is to be meant in the black hole case as 
the inverse black hole temperature. Still, we point out that the following analysis holds true for 
fermions in a generic thermal state with inverse temperature p. 

We refer both to and to |Q , and first we consider the Bogoliubov transformation for the "diag- 
onal" case (we purposefully choose a notation which allows a straightforward comparison with ]32| ; 
see also Q): 

ar* = Mi ar + '^i (bD^ (2.11) 

br' = br ~ {an\ (2.12) 
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where / is a collective index specifying states in the Hilbert space, and where the usual CCR rules 
for fermions lead to 

|^/p + |i.,|2 = l. (2.13) 
We are interested in the following operator: 
Ar°"*(a) : ^ (ap*)tap* 

= \^ii\\arVar + ^^v^rar + ^^*KoHb^y + h\' (i - m^br) . (2.14) 

For our aims, it works equally well the operator which allows to detect the net effect of the instability 
(cf. H) 

Nr\a) := (aP*)1'aP* - (ar)^ar- (2-15) 

We also introduce thermal state operators, according to the standard constructions in thermofield 
dynamics ||2^, |3^, |3^, both for initial and final states. We omit "in" and "out" labels in this case, 
for simplicity of notation, and introduce the thermal state |0(/3) > and thermal state annihilation 
operators ai{f3),di{P),bi{(3),bi{(3), which arc such that 

aimO{(3) &i{f3)\0{f]) >= bi{l3)\0{l3) >= bi{(3)\0{p) >= 0. (2.16) 

We are mainly interested in the following relations (see also p6t): 

ai = s+aiiP)+c+~4{f3), (2.17) 
bi = s-bi{p)+c-bl{P), (2.18) 



with 



c+ • = ^ (2 19) 



exp[l/3(a; - ^+)] 
VI + exp[(3{u; - <^+)] ^ 



(2.20) 



and the analogous ones for 6-operators (which correspond to operators for antiparticles, i.e. for 
negative frequency states; cf. p9[|): 

^ (2.21) 



^1 + exp[P{\uj\ + if-)] 

(2.22) 



exp[i/3(|^| + ^-)] 



y/1 + exp[(3{\uj\ + ^-)] 

where ip^, tp^ are chemical potentials for particles and antiparticles respectively |]36| . 
In the previous formula and in the following ones, for simplicity of notation we make a "liberal" 
use of indexes for quantum numbers, as far as unambiguous formulas arise. By making use of the 
above relations between particle operators and thermal state creation and annihilation operators, 
we easily find for particles 

< iVr* - {ctf K?) = Wil'l (^tanh[i/3(c.- ^+)] +tanh[i/3(|..| +y.-)]) (2.23) 

this result in the limit as (p+ , (/?^ -^0 agrees with the result displayed in |Q . Note that this means 
that 

< Nr >P= {ctf + - {ctf - {cif) ^ < Nout (2.24) 
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where the former term is the expected mean number of fermions in thermal equilibrium at the given 
temperature and the latter term is the net effect associated with the pair production induced by 
the presence of an electrostatic field. In our notation for the the Nariai case, we shall get 

< Nr' >/3= |Tfc(c^)pi (^tanh[i/3(c. - ^+)] + tanh[i/3(|c.| + ^-)]) (2.25) 

where (p~^ is assumed for definiteness to be the chemical potential for particles in the case of a 
positively charged black hole, ip~ = ip^ , particles are electrons (charge — e) and 

ip+ ^ -e{Ao\rr ~ Aa\o) ^ -2eQj. (2.26) 

Formally, an analogous experession holds for the ultracold I case, which is nevertheless pathological 
(cf. sect. 4^). Furthermore, in analogy with Page's analysis in one can expect that charged 
particles are efficiently emitted with thermal spectrum only for small black hole masses, and then 
that above the corresponding threshold charged particles are emitted only because of the electro- 
dynamic instability (i.e. as if they were emitted in vacuum). 

Moreover, we have to take into account that both in the ultracold I and in the Nariai case we work 
with dimensionless (rescaled) variables and then we get (3 = 27r. Note also that 

< m|A^r*|0 in >= \Tk{uj)\'^. (2.27) 
It is also interesting to evaluate the following quantity: 

(AA^;)' iNr*f >/3 -< >\- (2.28) 
As to the operator (Ar°"*)2, being (iV°"*)2 = one finds 

(AiVO' = < K"* >/3 (1- < K"* >/3) 

= (c+)2(l - (c+)2) + \vi\\\ ~ \vA^) - \vA\\ - \vi\^)Uf + (cD'] 
+ \vA%c^f + (cD'Kl - {c^f - {c- f) - A^i\\ctf{\ - (c+)2 - {c- f). (2.29) 
In a stable equilibrium situation, i.e. in our case in absence of electrostatic charge, one would obtain 

(AiV0^ = (c+)2(l-(c+)2), (2.30) 

which corresponds to the first contribution displayed above. Moreover, it is interesting to point out 
that 

< m|(iV°"*)2|0 in > -(< m|iV°"*|0 in > f = |i^,|2(l - (2.31) 

which amounts to the temperature-independent contribution in the previous formula, where also 
a third contribution, which "mixes" the pair creation effect to the thermality of the background, 
occurs. 



3. Ultracold II case 

The ultracold II metric is obtained from the Reissner-Nordstrom-de Sitter one in the limit of 
coincidence of the Cauchy horizon, of the black hole event horizon and of the cosmological event 
horizon: r_ = r_|_ = Tc. See |l5|. In particular, the metric we are interested in is 

ds^ = -dt^ + dx^ + ^{d9^ + sm^{9)d(t)^), (3.1) 

with x e M and t e K.. Then the {t, a::)-part of the metric is a 2D Minkowski space, to which a 
spherical part is warped with a constant warping factor. One gets r§3 = — sin(0) cos(0), Ffs — 
cot{6). The electromagnetic field strength is F = —^/Xdt A dx, and we can choose Aq = VAx 
and Aj = 0, j = 1,2,3. It is also useful to define E := \/A, which represents the intensity of the 
electrostatic field on the given manifold. We note that it is uniform, and then one expects naively 
to retrieve at least some features of Schwinger's result. 
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3.1 The transmission coefficient approach 



Let us consider the reduced Hamiltonian which can be obtained by variable separation as in |12| : 
starting from the full Dirac equation {p — fi)"^ = 0, variable separation leads to the following 
reduced Hamiltonian: 

-e^/Ax - dx + V2Ak ' 



hi, 



(3.2) 



fc = ±(j + i)GZ — {0} is the angular part eigenvalue, /i and e are the fermion mass and charge 
respectively. Then one obtains the following equation for <i> = e"*'^^^' = ' 



-{eEx + Lj)l2 + i<J2dx + V2Afccr 



$ = 0, 



(3.3) 



where <Ji are Pauli matrices and I2 is the 2x2 identity matrix. Let us take the unitary transformation 
$ = e"*^'^!^ so that ai ^ e~ 



*i°'^(Tie'4°'^ and in particular 

(12,0-1, (72, 0-3) ^ (l2,0-i,(T3,-cr2)- 



Then one obtains 

~{eEx + uj)l2 + i<^3dx + V2Afco'i + ^(72 
which amounts to the following couple of differential equations 



Then we can get 



V2Afc - in 

and the following equation for ^1 is obtained: 



-{eEx + uj)^i + idx^i + (V2Afc - ifi)^2 = 0, 
(%/2Afc + - {eEx + uj)^2 - idx^2 = 0. 



1 



[{eEx + - idxii] , 



+ [{eEx + Lo f + ieE - fil] ^ = 0, 



where /x^. :— fi^ + 2AA:^ is the effective mass corrected by K-K modes. We can define 



— {eEx + u), 



and then we obtain the following equation 



1 2 
4^ 



j4_ 

2eE 



(3.4) 
(3.5) 

(3.6) 
(3.7) 

(3.8) 
(3.9) 
(3.10) 



whose solutions are parabolic cylinder functions [Q. The calculation is completely analogous to 
the one performed in [p| , and as in S] one can easily show that the transmission coefficient satisfies 



\Tk{c^)\' 



exp 



eE 



(3.11) 



The latter expression corresponds to the mean number per unit time and unit volume of created 
pairs and coincides with the WKB approximation for the same coefficient [|l^ . This means that the 
WKB approximation is actually exact for the given case. This result is expected, being our case 
easily realized as a completely analogous to the standard case except for the compact character of 
the 2D transverse space. By adopting the strategy in 1^, or also the one in for determining the 
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factor preceding the exponential terms, one obtains the following imaginary part of the effective 
action: 



1 eES 



l^W = --^ ^ g(fc)log(l-exp(-^l), (3.12) 
where g{k) = 2{2l + 1) is the degeneracy factor. 

A comparison with the standard flat space-time case ^ shows that the structure of the ultracold 
II manifold, which is a product of a 2D flat-spacetimc times a 2D sphere, yields to the substitution 
of the (Gaussian) integral over the transverse dimensions with an infinite sum over fc, still to be 
valued. We recall that in the flat space-time case, one finds for the 4D density of the imaginary 
part of the effective action 

- = 2^E;iI-p(^-^nj, (3.13) 

n— 1 ^ 

which, by avoiding to perform the integration on the transverse variables ^j^, becomes 

n— 1 ^ 

' a p±^ log 1 — exp ' 



which is of course strictly analogous to (3.12). 
3.2 The C-function approach 

We tackle the problem also by the zeta function method. This will be useful later to compare with 
the more difhcult Nariai case. The Euclidean formulation is used. For convenience, we recall that 
the spectral zeta function associated to an operator H having eigenvalues A„ with degeneration c?„, 
is defined by 

that is substantially the Mellin transform of the kernel^ Kh{x) = Tre^^^. The point is that 
— logdetiJ = ^C-ff(O) defines the Euclidean effective action. After turning back to the Lorentzian 
signature, the instability is measured by the imaginary part of the effective action, which must thus 
compared with the vacuum persistence computed with the previous method. To this aim, let us 
consider the Dirac Euclidean equation for the ultracold II background. The spectrum of the Dirac 
operator, as well-known, is neither positive definite nor semi-bounded. To overcome this problem, 
one relates the zeta function of the Dirac operator^ to the zeta function of its square, and for the 
sake of completeness we recall some key points. To get the right definition, it is convenient to start 
with an heuristic reasoning. Let 5* an cigcnfunction of the operator _p — fi: 

{p - At)* = A*. 

Then we have (A± + /i) = ^^Jlp^ and thus we can formally write 

1, , r^.. 1, ... /^.. 1 



log(det(^ -i^)) = \ log(dct(-^ + v/|>2)) + ^ log(det(-A* - V^^)) = L iog(det(A.2 

The factor 1/2 arises from the double degeneration of each eigenvalue ^ . Thus, it is convenient to 
define 

-log(det(^- m)) = ^C2_^.(0), (3.14) 



■^We assumed here that the spectrum is strictly positive, but in general this can be weakened by complex analyticity 
techniques. 



'if ijf) - ti)'^± = \±^± then, for example, + = \+^±. 
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and then for the Euclidean effective action we get 

W = \c'^._p.{Q). (3.15) 
The simplest way to proceed is to compute the eigenvalues of —p^, and next to add the mass square 

We fully exploit the K-K reduction in order to perform the ^-function calculation. We first 
note that for the 4D Dirac operator we have 

lp=:$+r, (3.16) 

where f depends only on variables for the first 2D factor of the metric and analogously f depends 
only on spherical variables of the last 2-sphere factor. When one considers — -D^, it is easy to realize 
that one obtains 

and then the eigenvalue of —D^ is the sum of the eigenvalue uP' of —E"^ and of the eigenvalue 
&^fc^ of {b is a constant related to the radius of the 2-sphcrc factor and k is the usual eigenvalue 
for the angular momentum operator K, which is such that —F^ = b^K^): 

X^ = w^ + b^e. (3.18) 

In the ultracold cases, one finds 6^ = 2A; in the Nariai case, one has b^ = B. Eigcnfunctions for —D^ 
(and then also for —D^ + fi^) are tensor products of eigcnfunctions of —E^ and of eigcnfunctions of 
—F^. As a consequence, degeneracy can be read from the aforementioned tensor product structure. 

In what follows, we define -jfj,, /j. = 0,1, 2, 3, as the Euclidean version of the usual gamma ma- 
trices which are such that {7/:j,7/n} = ^S^^v 

In the present case, we obtain the following axpression for E: 

$ = lo{dt - ieEx) + ^A, (3.19) 

and then 

E"^ = {dt - ieExf + dl + ieE joli ■ (3.20) 

The first two terms are meant to be multiplied by 4D identity I4. Being iei?7o7i = eEa2 (S> I2, and 
I4 = I2 (8) II2, the eigenvalue equation 

-E'^ij = w^ip (3.21) 
is equivalent to the following reduced 2D equation 

[-{dt - ieExfh - d% - eEa2] C = w^^; (3.22) 

a further unitary rotation such that o"2 '-^ cts carries the problem in a "diagonal" form 

{-dl - {dt - ieExf T eE)ri± = w'^r]±, (3.23) 

where t]± are components of the 2D vector eigenfunction of the above eigenvalue problem. Variable 
separation 

V± {t,x) = exp(-iw<)C± {x) (3.24) 
leads to the following differential equations: 

dlC± + [±eE - (w + eExf + w^] = q- (3.25) 
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defining the new variable z 



we get 



C±W = o, 



(3.26) 



which is easily realized to be the equation one obtains for the standard simple harmonic oscillator. 
The quantization conditions are 



2n + l = ±l-— , 
eh 



which lead to 



(3.27) 



and then one can write 



— 2neE, 
w'i ^ 2{n + l)eE, 



w — 2neE, 



(3.28) 
(3.29) 

(3.30) 



by taking into account that n = has degeneracy which is one half the degeneracy of n > (cf. 
also the 2D case in [pof). 

As a consequence, the eigenvalues for —D^+^^ are 



= 2neE + fil, 



(3.31) 



where fif. = 2Afc^ + Notice that the eigenvalues do not depend on w. An overall degeneracy 
factor d must be determined. This can be obtained by comparing the behavior of the kernel K(x) 
in a; = to the universal coefficients provided by the heat kernel theorems. We have 



d [ 2 ^ exp{-2neEy) - 1 



n=0 



The part in square bracket is the heat kernel for the 2D operator —E'^, and in the limit as y ^ 
it holds [• • •] 2d(2eEy)-^. As it must be equal to 25(47ry)~i (see [|o|, p. 368), where S is the 
volume of the 2D space, we get d — . We now define 



C(s) =:^.g(A:)a(s), 



(3.32) 



with 



eES 



/ dtt'~^y2exp{-2neEt-fj:it)~-— dtt''^eM~At) 
Jo „ n ^ ^ y'^i J a 



2-K 



n=0 

(2e£;)-Xff(s,— ) - ^ 



(3.33) 



By putting eE i— > ieE wc finally obtain 

1 



i-k(o) = f^ 

A Ztt 

eES 



2-K 



log(2ei?) + ^ - ^ log(2^) + Rclogr (^-^^ j + - log(/.^) 



(3.34) 



The volume factor appears here because of integration over the whole spacetime is included. This 
result agrees with the previous one. 



*We introduce a factor i in view of (3.15) 
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4. The ultracold I case 



A second extremal limit of the Nariai background is given by the type I ultracold solution when 
r_ = r+ = Tc. The metric is jl5|l 

^ _;^2^^2 ^ ^^2 _^ -L.(d6»2 + sin2(6')d<?!)2), (4.1) 

with X G (OjOo) and '0 € M, and the electromagnetic field strength is = ^ ^ifj. The 

spacetime presents the structure of a 2D Rindler manifold times a two dimensional sphere (with 
a constant warping factor). One gets rQ;^ = x'^oo — X^^lz = ^ sin(6') cos(0), = cot(6'). We 
can choose = -^X^ ^-nd Aj = 0, j = 1, 2, 3 as potential. The situation is now a little bit more 
tricky, but we still are able to compare the transmission coefficient approach with the zeta function 
method. We expect to find out results which are analogous to some extent to the ones obtained in 
1^ for the case of a charged scalar field in a 2D Rindler spacetime with external electrostatic field. 

4.1 The transmission coefficient approach 

In we have obtained that variable separation allows to obtain the following reduced Hamiltonian: 



hk = 



-^X'-MX xdx + V2Akx 
-xd^ + V2Akx -^X^+I^X 



(4.2) 



Again /i, e are the mass and the charge of the fermion. We also introduce E a/A- 
Using the coordinate t ~ x^/2 the Dirac equation takes the (Hamiltonian) form p 



\{eEt + cj)/2 + [i\J2tG'i - i1ta2dt - 2V2Afccri]C = 0. 
Rotating to ^ by the a\ Pauli matrix as before: C ~ e^*T°'i^ , we get 



[2t5t + iieEt + w)]^i - (/i + iV2Ak)V2t^2 = 0, 
[2tdt - i{eEt + w)]^ - (m - iV2Ak)V2i^i = 0. 

This is equivalent to the system 

{^l + iV2Ak)V2t^2 = [itdt + i{eEt + cj)]Ci (4.3) 
= tCi + - ^ [2iMfe - (eEt + i{eEt - LoM^ . (4.4) 

Set £,i{t) = t^i^ e~^'^^* F[t). Then, F must satisfy the confiuent hypergeometric differential equa- 
tion ^ 

= tF" + {^-iuj- ieEt)F' - ^F, 
which has general solution in terms of Kummer functions $(a; c; z): 

Using the Kummer relation 

$(a; c; z) = e^$(c - a; c; -0) 



we get 

^i(t) = ae-i'^^H-i'^^i-^; I - iuj; leEt) + /3e5'=-^4^+5'^$(l - ^ + zw; -ze£;i). (4.5) 

2ieE 2 2ieE 2 
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To compute ^2 we use (4.3) and 



$'(a; c; z) = -$(a + 1; c + 1; z). 
c 



We get 



leEa 



^e-^'^-Bt^-f j4_ ^^Et) 

^ + i^Ek ieE{l - 2iuj) ^2ieE '2 ^ 

^l + i^/2Ek^eE{i + 2^uJ) ^ 2ie£:' 2 ^ 

2t{fi + iV2Ek) 2ieE'2 ' ' 



The asymptotic behaviors of these sohitions are, for t w 

V2 



^1 «ae-^"-^*r5'^, 



1, 



whereas for t 



-00 



eEt 



a{-ieE)-A £(2 + p{-ieEy 



(4.6) 

(4.7) 
(4.8) 



r(i 



2ieE > 



(4.9) 



^2 « es'^-^'ii^+A. ^ 



/i + iV2Ek 



r(5^) r(i+z..+ ^). 



(4.10) 



It is useful to introduce a new variable x :— log(x). Using group velocity, (4.9) and ( 4.10 ) represent 
ingoing and outgoing particles respectively at x ~ 00. Similarly, (4.7) and (4.8) are the outgoing 
and ingoing particle respectively aX x k —00. In this situation, we can compute the transmission 



and reflection coefficients. Rewriting (4.7) and ( |4.8| ) as 



and (U) and ( [4.10D as 



with (t){x 
get 



ie£;e2^, then we must put C so that R ^ A/ B and T = D / B. We 



cosh[7r(a;- ^)] 



e 2eE 



cosh(7rijj) 



cosh[7r(u; - ^)] 



(4.11) 
(4.12) 



As a check we note that |Tfc(w)p + |i?fc(cj)p = 1, and |Tfe(cj)p gives the mean number of created 
pairs per unit time and unit volume. 

A comparison with the result obtained in the WKB approximation in |^ is also in order. We 
recall that WKB gives for |Tp the same result as in the ultracold II case. Then, the dependence 
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on w is missing. We can relate this result with the exact one as follows: one considers eE fixed 
and a; ^ — oo. It is easily seen that the above exact result converges to the one of the WKB 
approximation in this limit. 



We now calculate the imaginary part of the effective action. We write again 



ImT4^ ■.= yg{k)Wk. 



k 



(4.13) 



First of all, we note that there is an important difference with respect to the type II case. Here, 
the level-crossing region, i.e. the region where particle and antiparticle states overlap (cf. e.g. [^), 
is no more the whole energy range but, assuming for definiteness eE > 0, is determined by a; < 0. 
As pair production is expected to happen only in the level-crossing region for the case eE > 0, we 
must calculate 

Wk ^ -lY.\og{l - \n{uj)n (4.14) 



for uj < 0. We have 
easy to show that 



^- j duj (cf. ^), where T stays for a finite time interval, and it is 



log(l- |r,n = log 1-exp 



eE 



log 1 + exp 2iTuj 



eE 



We have to evaluate the integral 



duj log 1 + exp 2ttuj 



eE I 



and then we get 



Wk 



1 r 



duo ) log 1 — exp 



eE 



2ti 



Li2 



exp 



eE 



(4.15) 



(4.16) 



(4.17) 



The factor ^ duj^ amounts to a degeneracy factor which we can evaluate following Q]. As 

explained in the introduction the geometry of the ultracold I manifold is of the form of a Cartesian 
product M^'^ X 5^, where the sphere could be considered an internal space. Thus, this situation 
can be treated in the same way as in ]44| where a 1 -I- 1 dimensional space time is taken into account. 
There, it is considered the Klein-Gordon equation for a complex scalar field in the presence of an 
E field, but the same reasoning also works for a Dirac field. In this way we obtain a Schrodinger 
equation for a particle in an upside down oscillator potential and we can compute the degeneracy 
factor counting the number of modes whose turning points lie within < "0 < T and < % < L, 
where T and L are the sizes of the space time box over which E is nonvanishing. We put TL — S, 
and then we obtain for the degeneracy factor the value of eES /2Tr, exactly the same as we will 



obtain by using the C- function approach. We underline that ( 4.17 ) is in strict analogy with the 
results obtained in |^5| for the scalar case in a 1-t-l Minkowski spacetime, and the terms appearing 



in (4.17) allow an analogous interpretation: the first contribution is leading and is a bulk one, 
proportional to the spacetime volume of the 2D Rindler part of the manifold, the latter one is a 
surface contribution. Cf. also [M. 



4.2 The ^-function approach 

The present case is also analyzed by means of C- function techniques, which confirm the results 
obtained in the former approach. 
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For the Dirac operator on the {tp, x)-part of the manifold one obtains 

¥ = ^70 (d^ - ieE^^ + 71 + l^y, (4.18) 
A Liouville unitary transformation 

(5*)(V,x):=x/X*(V',x) (4.19) 
(i.e. ^(V',x) = :^^(V')X))) leads to the unitarily related operator (still called^) 

¥ = ^70 (^v, - ieE^^ + 71 (4.20) 

and then ^ 

E' = ^ - y) +91+ 7071 {^eE^ + ^d^) (4.21) 



Substituting t = performing variable separation and keeping into account a unitary "rotation" 
transformation which is completely analogous to the one performed in the ultracold II case, one 
obtains for the eigenvalue equation of —E^ the following couple of ordinary differential equations: 

td^r]±{t) + ^dtT]±{t) - ^ [(w + eEtf ± {to - eEt) - 2tw''] rj±{t) = 0. (4.22) 



Let us introduce g±{t) through 



ri±{t) = exp ( ±^eEt]t^--'^g±{t); (4.23) 



then we obtain the following confluent hypergeometric equations: 

t 



i ± (w + eEt) 



dt 2 



We are looking for solutions ( '^^ ) which belong to the Hilbert space L^[(0. 00). inherited 

\V- J 

from the original 4D space (the measure keeps into account the above Liouville transformation). 
We need to distinguish two regions for the energy ui. 
For a; < 0, we find 

-w^ 1 



with the quantization condition 



5(i)-=^(^;2-^;e^^*)' (4-25) 



(4.26) 



2eE 
which gives 

w'i = 2eEn. (4.27) 
Moreover we find ^ 

g{t)+ = eM-^Et) $(1 - I^; ^ - w; eEt), (4.28) 



with the quantization condition 
which gives 



l-^ = -n, (4.29) 



wl = 2eE{n + 1). (4.30) 
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Moreover, as in the ultracold II case, we can re-label the eigenvalues in such a way that the 
eigenvalues of —p^ + p? in this region are 



(4.31) 



where again /i^ = /i^ + 2Kk^ . Also in this case the degeneracy of the n = case is one half the 
degeneracy of the n > cases, in full analogy with the ultracold II case. For the region w > 0, 



solutions ^ ^ G -^^^[(0, oo),yT correspond to 



and to 



.{t) = t5+-$(_ + ^ _ _; _ + ^; eEt), 



1 I 
g+{t) = exp{-eEt)^{- + w - — ; - + w; eEt), 



which are both associated with the quantization condition 



(4.32) 
(4.33) 



1 



2eE 



-n. 



The eigenvalues of —p"^ + jj^ are for a; > 



= (2n + l)eE + 2eEuj + (4- 



(4.34) 



(4.35) 



Note that the explicit dependence on iv appears only for the second region. For w < 0, i.e. in the 
level-crossing region, we have independence on iv of the integrand. The degeneracy factor can be 
computed as for the ultracold II case giving the same factor. As to the term depending on w, with 
w > 0, we recall that J2uj '"^ ^ / ^"-^ holds. Then we obtain 



2^*=^ ^ 27r 



\2eE)-CHi^,s) 



1 1 
^1^ 



T 

2^ 



{2eE)- 



^C.(^(l + §),.-l).(4.36) 



Note that the first term is the same as for the ultracold II case. Differentiating the second term 
with respect to s in s = and going back to the Lorentzian signature we obtain the contributions 



T 

2^ 



(log(2ieE) - 1)Ch(^[1 + + ^log(MD - C'Hilii + 



To compute the imaginary part of this expression, we first note that 

CH(a, -1) = -^S2(a) = -^(a^ - a + ^), 
(^2(0) stays for the second Bernoulli polynomial) so that if a is real 



(4.37) 



ImCff(-+ia,-l) = 0. 



Then, the first term gives the contribution 



1 



2\eE J "^6 



To compute the second term, let us start from the obvious identity 



(4.38) 



(4.39) 



(4.40) 
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Deriving in z = — 1 and choosing a = ix we get 



Cni^l + \: -1) = ^Ch(«; -1) - CH(*f ; -1) + \ log2Cff (zx; -1). 



Using l-mC,H{ix; —1) = | and formula (20) in G2| we get 



^1 , X 1 , 1 



_2 1 °o p-2nTr3: 

-Y-2("") 

n— 1 



E 



71^ 



2E- 

n=l 



1 

8^ 



2 1 ^ 



Using the expression 



E(-ir 



= Li2(-e" 



and adding ( 4.39| ) and the derivative of the first term in ( 4.36 ), we finally get 

1 



1 ^ eES 
-ImCUO)^^ 



'2 ( 1 - exp 



Ll2 

27r 47r 



exp 



(4.41) 



(4.42) 



(4.43) 



(4.44) 



This result coincides with (4.17). It is remarkable that the transmission coefficient approach involves 
only an integral over the level-crossing region, whereas the (^-function calculation requires a control 
over the whole spectrum, and is quite more tricky, but it provides also more information, being a 
priori the complete 1-loop effective action made available by the latter approach. 

4.3 Instability of the thermal state 

According to the general discussion of section 0, we get 



e 2eB cosh(7rijj) 1 , , . , , , , , -,^\ 

^ -- (tanh[7r(tj -</?+)]+ tanh[7r(|cj| + </j )]) 



l + exp[27r(tj-(^+)] cosh[7r(tj - ^)] 2 

(4.45) 

where the latter term is the net effect associated with the pair production induced by the presence 
of an electrostatic field. Notice however that the potential ip^ = (p~ is ill-defined, being infinite 
unless a spatial cut-off is introduced at x = Xo > 0. 



5. Nariai case 

We now consider the more general case, that is the electrically charged Nariai solution. The manifold 
is described by the metric jl^, |lj, |l^ 

ds^ = ^(- sin2(x)dV^ + dx^) + ^{dO^ + sm^{e)d<j)^) (5.1) 
A. B 



with e R, X G (0, tt), and the constants B = ^ |^1 - y 1 - 12^ j , A = - B aie such that 

^ < 1, and :— j^. The black hole horizon occurs at x = This manifold differs from the 
ultracold cases because it has finite spatial section. In the Euclidean version, it corresponds to two 
spheres characterized by different radii. One finds the following non-vanishing Christoffel symbols 
ro;^ = cot(x),rJo = sin(x) cos(x), r§3 = — sin(0) cos(0), — coi{9). For the gauge potential we 
can choose Ai = —Q^ cos(x)^i'- 

As the situation is more difficult, in this section we will provide a more detailed exposition. 
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5.1 The transmission coefficient approach 

By variable separation, as performed in | p^ , wc obtain tlic following reduced Hamiltonian (where 
/z, e are the fermion mass and charge, as before): 



eQf cos(x) - ^ sin(x) sin(x)5x + sin(x)fc 
- sin(x)5x + sin(x)fc cQf cos(x) + ^ sin(x) 



We introduce, with the aim of simplifying the notation, the following definition: 



(5.2) 



(5.3) 



which corresponds to ^ times the maximum value for the intensity of the electrostatic field. More- 
over, we assume for definiteness eE > 0, and also we adopt the following redefinitions: 



fj. ^ fj., 
k k. 



(5.4) 



Then the Dirac equation in the Hamiltonian form is 



eE cos sin X sin xdx + k sin x 
— sin xdx + k sin x eE cos x + sin x 



Using the coordinate t — — cos x we can write it in the form 

{eEt + uj)l2 + Mv/T^f^CTs - i{l - t^)(T2dt - k^Jl - f^ui 

where 5* — 



^ ~ (^"2) ■ tSike the unitary transformation 5* = e 'I'^i^; then 



From ( |5.q ) we find 



[(1 - t^)dt + i{eEt + u;)]^i - (^ + z/c) VT~^6 = 0, 
[(1 - t^)dt - t{eEt + uj)]^2 - (a* - tk)VT~t^^i = 0. 



.N^ / n ^, eEt + uj ^ 



which inserted into (5.7) gives 



(1 - - te; - (m' + k^- »ei?)Ci 



l-t2 



[(e£'i + cj)^ + iieEt + tj)<]^i = 0. 



(5.5) 

(5.6) 
(5.7) 

(5.8) 
(5.9) 



Looking at the behavior of this equation at the singular points t = ±1 we see that it is convenient 
to define a function C such that 



a(t)-(i-t)"(i + t)''c((i + o/2), 



(5.10) 



\ i \ i 

ae {-(l-ieS-iw),-(eS + w)}, /3 G {-(1 - ie£; + iw), -(e£; - w)}. (5.11) 

Setting z = (1 + i)/2 and choosing a — \(eE + uj) and /3 = §(ei? — w), wc find 



z(l - 0)C" + [2/3 + - - z(l + 2ieE)\C: - [fj.^ + = 0, 



(5.12) 
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which is an hypergeometric differential equation. Note that the solution regular in 2: = of this 



equation provides a solution of (5_^) regular in < = —1. Moreover, ( p. 9]) is invariant under the 
combination of complex conjugation with the transformation (t, E) (— t, —E). Thus, we are able 
to get the solution regular m t = I applying this transformation to the solution regular in t = — 1. 
This gives 



+ C2 



1 - t 



1-t 



1+t 



1+t 



2Fi{ieE + iVA, ieE - iVA; ^ + i{eE - w); ^-^) 
2Fi{ieE + i%/A, ieE ~ i^A; ^ + i{eE + uj); i— ^), (5.13) 



where A := /i^ + fc^ + e^£'^, and 2F1 is the well-known Gauss hypergeometric function. 
Using (|5.8|) and the relations 



1 - t 



1 + t 



eEt + UJ 



1 + t 



= 



ab 



2F[{a,b;c;z) = — 2-Fi(a + 1, 5 + 1; c + 1; z), 
c 



we get 



2(p~ik) fl~t 



l + 2i{eE-u;) \ 2 
; + iy/A + l,iej 
2{fi~ik) fl-t 



1 + t 



(5.14) 



3 1 + t 

2Fi{ieE + iVA + 1, ieE - iVA + 1; 2 + ^(^^ ^ ) 



C2 



1 + 2i(e£; + cj) V 2 



2 ^2 2 ^2 



2Fi{ieE + iVA + 1, ieE - iVA + 1; 2 + ^(^^ + 
As in the previous cases, we can look at the asymptotic behaviors at infinities, in the coordinate 



X = log tan ^. For x — > —00 we get 



r(i + i{eE + w))r(i - i{eE - u;)) 



r(i +iuj- i%/A)r(i + + ^^/A) 



(5.15) 



6(2;) ~ -C2 

For X ^ +00: 
a (a;) 



2{n~ik) r{^ + i{eE + uj))T{^ + iieE^uj)) 
1 + 2i(e£; + cj) r(l + ieE - iy/A)T{l + ie£; + i^/A) ' 



-i{eE—Lu)x _j_ ^i{eE—Lu)x0^^x 



(5.16) 



C2 + Cl 



r(| + z(e£; - u;))r(i - ijeE + u;)) 
r(i -iuj- iVA)T{^ ~iuj + iy/A) 



£ (a;) ^ c 2(/i-»fc) V{l+i{eE-uj))T{\ + i{eE + ^)) ^,^,e+u.)x 



(5.18) 



To compute the transmission and reflection coefficient let us write the asymptotic expressions for ^ 
as 



£,i{x)- « [ci + C2ao] e''"-^-'^)^ + e-'('=^-")^0(e^) 
6(a;)" » -C2 /Joe"*'^-^-^)^ + e*(^-^-'^)^0(e^), 



(5.19) 
(5.20) 
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for X — > —oo, and for a; — > oo 

^i{x)+ « [C2 + ciai] e-^(^-^+")" + e*(^^+")"0(e-"), (5.21) 

By taking into account correctly group velocity and imposing that there is only incoming wave at 
X — — cx), one finds 

2 cosh[7r(e-B — w)] cosh[7r(ei? + w)] 

\-lk[UJ)\ - 



cosh[7r(v A — Lo)] cosh[7r(v A + uj)] 
cosh[27rei?] + cosh[27rijj] 

[o.Zo 1 



cosh[27rvA] + cosh[27raj] 

which has the expected property to satisfy |Tp < 1 and gives the mean number of created pairs for 
unit time and unit volume. Moreover, we get 



cosh[7r(\/ A — uj)] cosh[7r(v A + uj)] 
and the property \Tk{uj)\'^ + \Rk{uj)\'^ = 1 is verified. 

We recall that one has to reinstate in the above formulas for jTp the original values for /i and k 
{ci. eqn. Q). Then A = ^ + f fc^ + (eE)^. 

Comparing with p^ , the limit as eE oo leads to the WKB approximation (and to the limit 
|r/c('^)P 1 )• This is the actual approximation where the WKB approximation works well. We 
can notice that the above limit actually means Q-? ^ oo. Being 



' A 

B _ 1 

A ~ ^1 - 4AQ2 



(5.25) 



one can obtain the above limit as - — 



4A; 



In order to determine the imaginary part of the effective action, we refer again to ( 4.13 ). We 
do not perform the sum over k, and then we calculate 

Wk^-lY.^ogil-\n{u;)f). (5.26) 

UJ 

We sum only over the level-crossing region, because only there particle creation is expected to be 
present, and then only there an instability for the vacuum should occur. This region is —eE < to < 
eE so that we have to perform the following integral: 

, , cosh[27re£;] + coshpTTO)] ^ 

dwlog(l ^ -), (5.27) 

-eE cosh[27rv A] + cosh[27ra-'] 

where the dependence on k is implicit in A; the integral can be rewritten as follows: 

/ = 2e£; log(cosh[27r%/A] - cosh[27re£:]) - //, (5.28) 

where 

dwlog(cosh[27r\/A] + cosh[27rw]) = — / dy \og{p + cosh[y]) , (5.29) 

-eE 27r J-2TreE 
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with p := cosh[27r'\/A]- By taking into account that p + cosh[?/] = \{2p + exp(?;) + exp(— y)) = 
\ exp(y)(exp(- 2y) + 2 pexp(-y) + 1) = i exp(y)(exp(-y) + p - ^Jp"^ - l)(exp(-j/) + p + ^/p^ - 1) 
and that p — yp^-l = cxp(— 27rVA) and p + yp^-l = cxp(27rVA), one finds 

r.27re_E 



27r 



27re£; 



y — log 2 + log(exp(— y) exp(— 27rv A) + 1) + log(exp(— y) exp(27rv A) + 1) . 

(5.30) 



The following result is useful: 

j dy \og{exp{- y + S) + 1) =iy2 + y log(exp(-y + (5) + 1) - ylog(exp(y - <5) + 1) 



Li2(-exp(y - 6)) 



1 



y^ + yS - Li2(- exp(y - S)). 



(5.31) 



As a consequence, with simple manipulations, we get 
1 



// = 



27r 



-2{2TreE) log2 + Li2(- exph27r(VA + eE)]) - Li2(- exp[27r(VA + eE)]) 



Li2(-exp[27r(\/A - e£;)]) - Li2(- exp[-27r(\/A - e£;)]) 



(5.32) 



and then, by taking into account that i— > ^ J dco, we get 



Wk 



T T 1 

e£; log(2 cosh[27r\/A] - 2 cosh[27re£;]) 

27r 2tt An 



Li2(-exp[-27r(VA + e£;)]) (5.33) 



Li2(- exp[27r(\/A + e£;)]) - Li2(- exp[27r(\/A - eE)]) + Li2(- exp[-27r(%/A - eE)]) 



In the Nariai case, it is not so straightforward to distinguish between bulk and surface parts of 
the given Wk- By keeping into account that E (x -j and that the volume of the (Euclidean) 2D 
(-0, x)-part of the metric is one realizes that the first above contribution is leading and is a bulk 
one and the latter is again a surface contribution. 

5.2 The function approach 

This case requires more efforts in order to be tackled with ^—function techniques, in particular 
it needs techniques which have been developed only recently pc| ], and we sketch herein a more 
heuristic approach which still leads to the desired results. 



As to the operator on the (7/;,x)-part of the manifold, one obtains 



smx 



70 {d^ + ieE cos x) + v A71 ( 9^ + - cot x 



(5.34) 



After a Liouville unitary transformation '^{tp, x) — ^slnx '^^^' ■^^^ simplified form (again 

called |;) ""^ 



smx 



70 {d^ + ieE cos x) + V Ajid^ ; 



moreover, 
E^ 



. \ {d^ + ieE cos x)"^ + Adl + A7071 
sm X 



cos X 

2 — (^V) + '-^'^^ x) + i^E 

sin X 



Substituting t = — cos x one obtains 



E^ = 



{d^ - leEt)^ + A{1 - e)d1 - Atdt + Ajoli (- " *^^^) + *^-^) 



(5.35) 
(5.36) 

(5.37) 
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After variable separation, a reduction and a rotation as in the previous cases, we obtain for the 
eigenvalue problem of —E^ the following couple of differential equations: 

1 

(1 - t^)d'^r]± - tdtr]± ± eET]± + [-{cu + eEtf ± t{uj + eEt)] tj± + —tj± = 0. (5.38) 

We choose 

77+(t) - (l-t)^(l + t)^.g+(0, (5.39) 

77_(0 = (l-t)"^(l + i)"^5-(i)- (5.40) 

Then by choosing t — 2z — \ we obtain the following couple of hypergeometric equations: 

.(1 - ^)%|^ ^U-.^\- i2eE + l)z) ^ + ^,,(.) . 0, (5.41) 



with a+ = eE+,J^ + {eEf,b+ = eE - if + {eEY,c+ ^ eE ~ uj + \, and 

z{l - ^)^^ +{-eE + u: + - + {2eE - l)z\ + —g^{z) = 0, (5.42) 



with a_ = + (e£;)2,6_ = -eS - ^/^^ + (e£;)2, c_ = -e£; + cj + i. We are looking for 

solutions ^ ^ G i'^[(0, 1), ^(Iz^]^' difhcult to realize that this condition depends on cj, 

and three regions can be identified. 
One can show that 

g+{z) ^ 2Fi{a+,b+;c+;z) (5.43) 

with the quantization condition 

b+ = -n, (5.44) 

which leads to 

wl = A{eE + n)'^ - A{eE)'^ , (5.45) 

together with 

<?_(z) = z^-''- (1 - zY—'^''-+''-\Fi{1 - a_, 1 - 5_; 2 - c_; z) (5.46) 
with the quantization condition 

1 - a_ = -n, (5.47) 

which leads to 

= A{eE + n + 1)^ - A{eE)'^ , (5.48) 



correspond to solutions of ( 5.41 ), ( 5.42 ) respectively which allow to get 
in the interval —eE < uj < eE. We can re-label the eigenvalue and obtain for —D^ + ij? 

= A{eE + nf - A{eEf + ^? + Bk^ . (5.49) 

Also in this case, one has to take into account that degeneracy for n = is one half the degeneracy 
for n > 0. 

In the region oj > eE one can show that the couple 

.g+(z) = zi-=+2Fi(a+-c+ + l,6+-c+ + l;2-c+;z) (5.50) 
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and 

g^{z) = (1-2)^— («-+''-)2Fi(c_ -a_,c_ -6_;c_;z) (5.51) 
under the quantization conditions &+ — c+ + 1 = —n and c_ — a_ = — n, which both lead to 

= A(^uj+^+r?j - A{eEf (5.52) 

provides a solution which is in i^[(0, 1), ;j^^z^]^- Then the eigenvalue of —D"^ + fj,^ becomes 

1, 



X"^ =A{Lo + n+ -f - A{eEf + /x^ + Bk^ . 



(5.53) 

There is still the region w < —eE to be explored, where 

g+{z) = (1 - ^)'=+-(«++*+)2Fi(c+ - b+, c+ - a+; c+ - (a+ + 6+) + 1; 1 - z), (5.54) 

and 

g-{z) = z^'''^ 2Fi{l + b- - c_,l + a_ - c_;a_ + 6- + 1 - c_; 1 - 2) (5.55) 
with the quantization conditions c+ — a+ = —n and 1 + 6_ — c_ = — n, which both correspond to 



w'^ = A{-u)+^+n] -A{eE)'\ 



(5.56) 



satisfy the property to be eigenfunctions of the operator —E"^. The eigenvalues of — + /x^ are 

1 



= A{-oj + n+-f- A{eEf + /x^ + Bk'^. 

Note that, for |a;| > eE the eigenvalues of —D^ + /x^ can be written as follows: 

= A{\uj\ +n+^f+tx^ + Bk'^ - A{eEf, 

and then the integration for |u;| > eE is symmetric. 
For the heat kernel we get 



(5.57) 



(5.58) 



(5.59) 



with 



A({eE + nf + i^-{eE)'' 



Kk{s) = 2— {2 Vexp 
[ JeE ^0 

+ 2e£'^exp 

71=0 

where n\ = ^i'^ + Bk"^. Correspondingly, we obtain 



A[{^ + \+nf + ^-{eEf]s 



eEe^^{-fjils)) , (5.60) 



leE „ A" 



n + ^+uj 



1' + 4 - (e^)' 



+ (2e^) 



n A^ 



{n + eEf + ^ - {eEy 



1 1 
24" 



(5.61) 
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It is convenient to introduce the functions 



A' 



so that 



2eE ak{s 



--)-- — 



(Tfc(s; z)dz 



(5.62) 



(5.63) 



It is also useful to define a = ^ + eE and = ^ — {^EY . Using the Abel-Plana formula^ we get 




dx 



yl^ 



{x + a + zf + /52 



(5.64) 



{ix + a + z) + (3^ 



A' 



{-ix + a + z) + 



- 1 



To compute the effective action, we need to compute the derivative of (Jk{s\ z) with respect to s, in 

— . (5.65) 



s = 0. For the last term of (5.64) we get 

dx {a\\{ix + a + z)^ + - ln[(-ia; + a + z)^ + 



Note that this coincides with 



ds 



dx- 



1 



A'<^ {{ix + a + z + iPY {-ix + a + z + i^Y 



[ix + a + z — iPY {^ix + a + z — i^Y J e^'^^ — 1 



(5.66) 



where 7 is an arbitrary constant. The result does not depend on 7 and we could choose it equal 
to 0. However, it is convenient to choose 7 = ^- Applying Plana's formula to this integrals we see 
that 

d . , . d 



ds' 



=ocrk[s; z) 



ds 



.oCTfc(s; z). 



(5.67) 



where 



1 



^kis;z) =- 



2 A" \ {a + zf + f]^ 



+ 



^ 00 

7/2 X/ 



dx 



A" 



1 1 



{x + a + zy + ^2 
1 

(n + a + z + iPY {n -\- a + z — i(3Y 
1 1 



1 



2 A''/'^ \ {a + z + ipY (a + z ~ ipY 



As/2 



dx 



1 



1 



{x + a + z + iPY {x + a + z- if3Y 



(5.68) 



We now note that, thanks to our choice for 7, collecting the second and the last terms under the 
integral and deriving with respect to s in s = we obtain a vanishing term (see Lemma 1 in [po[ ). 
The same happens for the first and the fourth term so that only the third term contribute to the 
derivative. Our conclusion is that 



g-J.Mis;z)^-U^„ 



1 



As/2 



iCnia + z + if3,s) + Qnia + z — i(3, s)) 



(5.69) 



5 Abel-Plana formula is: Er=o /(") = 1/(0) + /o°° H^)'^'' + * Jo" -^'"t/ii"' ■ 
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Using d^el) we conclude that 



where 



with 



Thus 



d d ~ 

Tr\s=aCk[s) = -^\s=Q(^k[s), 
OS as 



2eE crfc(s;--) - --^ 



(Tfe(s; z)dz 



a-fe(s; z) 



{CH{a + z + if3,s) + (nia + z ~ if], s)). 



(5.70) 
(5.71) 
(5.72) 



T 
2^ 



2eE- 



Ai 



Cnia -^+i(3,s) + Cnia -^-if3, s) 



eE 



2^^^ [Cnia + ^/3, s - 1) + Cnia ~if3,s- 1)] 



A2 s -1 

Then we get (we momentarily omit the index k from some formulas below) 



(5.73) 



1 



T 



^C'(O) - ^^2(eS)2logA-e^logA + 2eSlog 



T{eE + if3)r{eE -i(3) 
2^ 



+ eSlog(^l) + (2 + log A) [Cff(a + iP, -1) + Cff(a - */3, -1)] 
- 2 [Cnia + 1(3, -1) + C^(a - -1)] 1 , 



(5.74) 



where we used the well known relations 



(,H{a,Q) = \-a, C^(a,0)=log^. 



(5.75) 



By going back to Lorentzian signature, through eE i— > leE, we get 

ilmC'(O) = ^ [-e£;iogA- e£;iog(A - (eS)2) - e£;iog f2cosh[27rVA] - 2cosh[27re£:' 
2 27r L V 

+ eE\og{^il) - 2Im [C^(a + ip, -1) + (^(a - -1)] 



(5.76) 



To compute the last two terms we can start from ( 4.42 ) and use the identity 

^2 



to obtain the relation 



and then 



-\\og\^z) - ^ = Li2(z) + Li2(l/Z), 



ImCH(*| + \-^)-^ [Li2(-e-^) - Li2(-e--^) 



ilmCfc(O) ==— (-e£:iog f 2 cosh[27r\/A] - 2 cosh[27re£;]' 
2 27r L V J 

\-U2{- exp[-27r(\/A + eE)]) + Li2(- exp[27r(VA + eE)]) 
-Li2(-exp[27r(VA- e£;)]) +Li2(-cxp[-27r(VA- eS)]) | . 



(5.77) 



(5.78) 



(5.79) 



Also in the Nariai case, the result coincides with the one obtained in the transmission coefficient 
approach. The calculation in the ^-function approach is much more difficult than the one in the 
transmission coefficient approach. Still, it furnishes the complete 1-loop effective action, and not 
simply its imaginary part. 
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5.3 Instability of the thermal state 

Also in this case, according to tlie general discussion of section |2[ we get 

1 



1 + exp[27r(tj - ip+)] 

coshpTreS] + cosh 27rtj] 1 / , , , i,, , r /i i m\ 

^ ^ -- (tanh[7r(w-(^+)] +tanh[7r(|cj| +</3")]) , (5. 



cosh[27r\/A] + coshpTTw] 2 
with (y9+ = 2eE ~ ip^ . Thermality of the state affects the pair production induced by the presence of 



an electrostatic field, which is associated with the second term in ( 5.80 ). We recall that in terms of 
physical (dimensionful) variables, by taking into account that = and that uophys — \fAw, 

in such a way that (3physOJphys = Zttw. 

6. Conclusions 

We have studied the spontaneous emission of charged Dirac particles by three special dS black hole 
solutions which share the relevant property to allow exact calculations and to present spherosym- 
metric metrics with a two-dimensional spherical part completely factorized (with at most a constant 
warping factor). As a consequence of the latter feature, the problem allows a reduction a la Kaluza- 
Klein to a two dimensional effective theory. This fact reflects itself in the common structure of 
our results in the ^-function approach as well as in the transmission coefficient approach, with the 
4D imaginary part of the effective action appearing as a sum over K-K modes of 2D terms. As to 
vacuum instability, our double check by the aforementioned approaches leads to identical results. 
Analogous results occur in the case of scalar fields on the same backgrounds |47|. Moreover, both 
the ultracold I case and the Nariai case make evident that calculations of the imaginary part of 
the effective action are easier in the transmission coefficient approach; still, it is to be pointed out 
that, a priori, the ^-function approach furnishes the complete 1-loop effective action (and then 



also vacuum polarization effects are taken into account |39|), and not only an evaluation of the 
vacuum instability (which is the only outcome of the transmission coefficient approach). It is also 
remarkable that the analysis of ultracold I and the Nariai cases points out that different regions in 
u) must be taken into account in the two approaches: in the transmission coefficient approach, only 
the level-crossing region is involved, instead in the C-function approach the whole spectrum must 
be included. 

As expected, the presence of an electrostatic field associated with the black hole charge induces a 
non-zero imaginary part of the effective action, i.e. to a vacuum instability which gives rise to the 
emission of charged pairs by the black hole, with a consequent discharge. 

The ultracold II case is to some extent trivial, because it substantially reduces to the instability of 
the vacuum in a 2D Minkowski space with an homogeneous electrostatic field. This fact appears 
clearly by comparing our results e.g. with the ones for the 2D case in Nevertheless, the results 
obtained are meaningful, as we deal with the instability of a black hole with zero temperature. 
The ultracold I case presents a less trivial structure, a non-homogeneous electrostatic field coupled 
with a Rindler-like horizon, and a non-zero temperature. Instability of the Boulware-like vacuum 
has been shown, and moreover also the instability of the Hartle-Hawking thermal state has been 
displayed; general calculations carried out in section |^ show that the transmission coefficient which 
signals vacuum instability still plays a relevant role for a thermal state at finite temperature, and 
Thermofield Dynamics is in agreement with the results in p^ . 

The Nariai case is the most interesting one and also the most difficult, because it represents a 
charged black hole characterized by a finite spatial section and with an electrostatic field which 
reach a maximum between the black hole horizon and the cosmological event horizon (where it 
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vanishes) , and then in this respect it is really different from Reissner- Nordstrom black holes (where 
the electrostatic field is mostly intense at the black hole horizon). Instability is again shown to 
occur. 

For the given manifolds, a comparison between the exact results for the transmission coefficient and 
the WKB approximation ones obtained in has also been performed, and, as expected, we have 
found that WKB is exact in the ultracold II case. 
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